clc 
clear all;
N=10;
m=4;
Va=randperm(N);
Vb=randperm(N);
Vab=[Va;Vb];
[newVa,newVb]=aPMX(Va,Vb,m);
newVab=[newVa;newVb];
disp(Vab);
disp(newVab);